Can the jamming transition be described using equilibrium statistical mechanics? 
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When materials such as foams or emulsions are compressed, they display solid behaviour above 
the so-called 'jamming' transition. Because compression is done out-of-equilibrium in the absence of 
thermal fluctuations, jamming appears as a new kind of a nonequilibrium phase transition. In this 
proceeding paper, we suggest that tools from equilibrium statistical mechanics can in fact be used 
to describe many specific features of the jamming transition. Our strategy is to introduce thermal 
fluctuations and use statistical mechanics to describe the complex phase behaviour of systems of 
soft repulsive particles, before sending temperature to zero at the end of the calculation. We show 
that currently available implementations of standard tools such as integral equations, mode-coupling 
theory, or replica calculations all break down at low temperature and large density, but we suggest 
that new analytical schemes can be developed to provide a fully microscopic, quantitative description 
of the jamming transition. 

PACS numbers: 05.20.-y, 64.70. Q-, 45.70.-n 
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From the point of view of statistical mechanics, the 
jamming transition observed by compressing random 
packings of soft repulsive particles in the absence of ther- 
mal fluctuations is an intriguing phenomenon [1]. It is a 
phase transition, in the sense that the mechanical re- 
sponse of the system changes abruptly at a critical den- 
sity [2]. It is also a critical phenomenon, since several 
scaling laws and diverging length scales have been de- 
scribed on both sides of the transition [3, 4]. Finally, it 
is a nonequilibrium phenomenon because the transition 
occurs in samples prepared out-of-equilibrium in the ab- 
sence of any relevant thermal fluctuations, and it is thus 
not possible to describe jamming without stating pre- 
cisely the protocol used to prepare the system [5]. Yet, 
quite remarkably, many features of the jamming transi- 
tion appear to be protocol-independent. 

Given the central role played by far from equilib- 
rium critical phenoncmona in statistical mechanics [6], 
jamming is thus a very active field of research and at- 
tracts the attention of the statistical mechanics commu- 
nity [1, 7]. In this conference paper, we mainly discuss 
the problem of the jamming transition from the point of 
view of statistical mechanics. We review, and expand in 
some places, recently published work on the subject, we 
discuss the philosophy and first results of our on-going 
effort to develop a new analytical approach to the prob- 
lem [8]. 

This paper is organized as follows. In Sec. I, we give a 
short account of the properties of the jamming transition, 
that theory should explain and reproduce. In Sec. II 
we briefly review the main theoretical approaches and 
explain the basic ideas behind our approach. In Sec. Ill 
we show how liquid state theory fares, and explain why 
it fails at low temperature. In Sec. IV, we treat the fluid- 



glass transition of harmonic spheres using both mode- 
coupling theory and replica calculations. We present our 
conclusions in Sec. V. 



JAMMING AS A NONEQUILIBRIUM 
PHASE TRANSITION 



In a pioneering study, Durian showed that jamming 
could fruitfully be studied in computer simulations us- 
ing simple models of soft repulsive spheres, and he intro- 
duced a model of harmonic spheres interacting through 
the simple pair potential [2] 



V(r) = e[l--) 0[1--), 



(1) 



where a represents the particle diameter, e is an energy 
scale, and 0(x) is the Heaviside function. Thus particles 
rcpcll each other harmonically when they overlap, but ig- 
nore each other otherwise. In Durian's original work, the 
potential was meant to describe the physics of wet foams, 
but this 'bubble' model is in fact so generic that it could 
equally be applied to the physics of soft colloids (such as 
dense microgels [9]), or emulsions [10]. Therefore, e can 
be interpreted as a parameter accounting for the elastic- 
ity of the soft particles, be they soft bubbles, colloids or 
droplets. In the absence of thermal fluctuations or exter- 
nal forcing, the unique control parameter for the phase 
behaviour of the model of harmonic spheres is the number 
density, p = N/V, for a system composed of N particles 
enclosed in a volume V. Equivalently, one can use the 
'packing fraction', tp = ira 3 p/6, although this name only 
strictly makes sense when particles do not overlap at low 
enough density. Indeed, because particles are soft, tp is 
not bounded. 
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FIG. 1: Packings of soft repulsive spheres at T = undergo a 
nonequilibrium jamming transition at a critical packing frac- 
tion tp c . The mechanical properties change from fluid to solid 
across tp c - Energy, pressure, and mechanical moduli increase 
continuously from zero above ip c , while the number of contact 
per particle jumps discontinously to a finite value at ip c . 



In Fig. 1 we summarize schematically the main fea- 
tures of the jamming transition that have been discov- 
ered, mostly through numerical observations [1]. Simula- 
tions have revealed the existence of a critical packing frac- 
tion, tp c , below which the packings have no overlap. In 
this phase, the energy, pressure and number of contacts 
between particles are zero. Above tp c , pressure, energy, 
number of contacts, bulk and shear moduli are non-zero. 
The linear (elastic) response of the system to deforma- 
tion changes abruptly at the transition: below cp c , small 
deformations can be imposed without any energy cost, 
while above tp c an arbitrarily small deformation causes 
an increase in energy. 

Interestingly, almost all these quantities increase con- 
tinuously from zero across (p c and vary algebraically with 
the distance to the transition, <p — ip c . A notable ex- 
ception is the number of contacts, which, as suggested 
by Fig. 1, jumps discontinuously to a finite value at (p c 
given by z c = 2d where d is the space dimensionality. 
The critical z c corresponds to the "isostatic" value, i.e. 
the minimal value for the system to be rigid [11]. At tp c 
the system is therefore marginally solid, with vanishing 
shear and bulk moduli. Moreover, the pair correlation 
g(r) of systems at ip c is very different from both amor- 
phous glasses or dense liquids, and present a number of 
singular behaviours, from a diverging contribution at the 
interparticle distance r = a [12, 13], to large-distance 
anomalies [14] reflected in a peculiar low-wavevector be- 
haviour of the isothermal compressibility [15]. 

Finally, the nonequilibrium nature of the transition is 
clear from the fact that thermal fluctuations are actually 
irrelevant for the features described above. This implies 
that 'crossing' the transition by 'compressing' the system 
is not a uniquely-defined procedure [16]. It is thus crucial 
to also specify how this is done in practice, in particu- 
lar how averages at each density are performed, that is, 
how different configurations at the same density are sam- 
pled. In an equilibrium setting, this is not necessary as 
configurations are sampled with their associated Boltz- 
mann weights. Two well-studied procedures to study the 
jamming transition are the following. First, rapid com- 



pressions of hard sphere systems simulated by molecular 
dynamics [17] produce jammed configurations in the limit 
of infinite pressure [18, 19]. Note that in this procedure 
thermal fluctuations play a role since thermal equilibrium 
can be reached at low enough density. Note also that the 
jammed phase cannot be accessed since particles cannot 
overlap. A second and very different procedure consists of 
studying systems of soft repulsive spheres, as in Eq. (1), 
directly at T = using energy minimization methods [3] . 
In both cases, sampling and averages are performed by 
repeating the compression or minimization protocol from 
a different set of initial conditions. A critical volume frac- 
tion tp c with the properties described above is generically 
found for any of these protocols. At present, there is nu- 
merical evidence that if the exact value of tp c is protocol 
dependent, the critical properties are not [20]. 



II. ON THE THEORY OF THE JAMMING 
TRANSITION 

At the theoretical level, there is at present no accepted 
framework to understand and account for all the fea- 
tures of the jamming transition that we have described, 
which mostly stem from very detailed numerical observa- 
tions [1] . Thus, if the physics and the main features of the 
transition are well described [7] , the theory is by compari- 
son less advanced. An elegant scaling theory, that identi- 
fies a divergent length scale and relates various observed 
scaling laws or critical exponents, was elaborated [4, 7]. 
Additionally, several distinct statistical frameworks were 
also developed [21-23], using in particular the idea that 
a statistical ensemble (the 'Edwards ensemble' [24] or a 
'force ensemble' [25]) distinct from the Gibbs ensemble 
must be introduced to study the jamming transition. A 
detailed account of these approaches is behind the scope 
of this short paper and we refer the reader to the recent 
reviews [1, 7] for further references. 

Here, we shall argue that the three main characteris- 
tics discussed above, namely the existence of a fluid-solid 
jamming phase transition, its associated critical proper- 
ties, and its nonequilibrium nature can in principle all 
be accurately computed within the standard framework 
of equilibrium statistical mechanics, starting from the 
sole knowledge of the interaction between the particles 
in Eq. (1). The general strategy we propose is to study 
the statistical mechanics of a system of harmonic spheres 
in the presence of thermal fluctuations. If this procedure 
is properly implemented, then we expect that a sharp 
phase transition with the above characteristics will occur 
at a critical volume fraction in the limit of T — > 0. 

The main conceptual difficulty to be faced has actually 
a physical origin, and deep consequences. In the relevant 
regime, the system is in fact characterized by the exis- 
tence of a large number of amorphous metastable states 
separated by large free energy barriers [5]. This means 
that the system gets naturally dynamically arrested in 
nonequilibrium glassy states, and that the free energy 
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landscape is quite complex. This implies first that the 
above mentioned dynamical protocols to probe the jam- 
ming transition are naturally affected by the glass tran- 
sition [5, 19]. Second, this means that statistical me- 
chanics treatments must accurately take the complexity 
of the phase space into account — or will necessarily fail. 
A theoretical framework to handle this complexity was 
developed in Rcf. [26] and applied in the context of par- 
ticle systems in Refs. [27-29]. It was also more recently 
implemented in Ref. [19] for hard spheres. 

The aim of the present work is to show that, although 
the approach is conceptually simple and direct, there 
remain purely technical difficulties, since previous im- 
plementations of these ideas are not very accurate in 
the limit of interest to describe the jamming transition, 
where temperature is low and density is large. In a sepa- 
rate work [8], we describe our ongoing efforts to circum- 
vene these technical difficulties to try and yield a predic- 
tive, quantitative, microscopic approach to the jamming 
transition. 



III. STATISTICAL MECHANICS OF 
HARMONIC SPHERES: FLUID PHASE 

A. Integral equations 

Our first step is to start investigating the behaviour 
of repulsive harmonic spheres in a regime where temper- 
ature is large. In the absence of attractive forces, the 
system is thus in a fluid phase, where correlations be- 
tween particles are small. 

In this regime, it is enough to treat the system using 
tools developed to study the statistical mechanics of the 
liquid state, such as integral equations [30]. These are 
typically closure relations yielding the pair correlation 
function of the fluid, 




where brackets denote a thermal average. 

In Rcf. [31] we studied the potential (1) using one par- 
ticular closure relation, the Hyper-Netted Chain (HNC) 
approximation. In this approximation, 

g(r) = exp[-pV(r) +g(r)-l- c(r)] , (3) 

where j3 = 1/T, c(r) is the direct correlation function 
defined through the Ornstcin-Zcrnikc equation: 

g(r) - 1 = c(r) + p J dr'cflr - r'\)[g(r') - 1]. (4) 

By numerical integration of Eq. (3), one can get g(r), and 
thus study the thermodynamic behaviour of the model 
for any state point (<p,T). 

Not surprisingly, no phase transition is detected by 
compressing the system at constant temperature within 
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FIG. 2: Nonmonotonic evolution of the pair correlation func- 
tion when volume fraction is increased at constant temper- 
ature in the fluid phase, as predicted using the hypernet- 
ted chain approximation [31], a simple integral equation fre- 
quently used in liquid state theory. Volume fractions are 
<p = 0.600, 0.663, 0.777, 0.900, and 1.00 (from right to left). 

such an approach, even when temperature is small. In 
fact, it can easily be shown that the ground state energy 
density, 

e g sO) = Hm e(T,<p) (5) 

remains zero at all p, since the energy density vanishes as 
e(T, tp) ~ T 3 / 2 , a result which is clearly inconsistent with 
the cartoon in Fig. 1. Simultaneously, the pair correlation 
function remains smooth and shows none of the jamming 
singularities described in the introduction in the range 
where the jamming transition is found numerically. As 
we argue below, this is not an artefact of the specific 
approximation employed, but is likely a generic feature 
of integral equations developed for liquids. 

Although a sharp jamming transition is not found, 
there is an interesting feature which emerges from the 
study of integral equations, as shown in Fig. 2. Upon 
compression at constant temperature, the evolution of 
the first and second peaks in the pair correlation function 
shows two distinct regimes, depending on the value of the 
volume fraction. For p < p*(T), the position of the first 
peak of g(r) shifts to smaller distances, reflecting the fact 
that particles get closer to each other. The height of the 
peak increases and structural order increases in the fluid. 
For p > p* , the position of the peak continues to shift to 
smaller r, but the height of the peak now decreases with 
p. Thus, the fluid becomes more disordered as density is 
increased, an 'anomalous' behaviour not seen in simple 
liquids. 

This density anomaly is actually well-known in the 
context of the physics of ultrasoft particles [32, 33], i.e. 
particles interacting with a pair potential that remains 
finite when particles fully overlap, as in Eq. (1) where 
V(r = 0) = e, but it is interesting to rephrase its phys- 
ical explanation in the context of jamming [31], where 
the crossover p* plays the role of a 'soft jamming' [34] 
or a 'thermal vestige' [9] of the jamming transition. It 
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FIG. 3: Numerical phase diagram for a binary mixture of har- 
monic spheres [36]. The symbols correspond to 'fluid' state 
points where thermal equilibrium can be reached using molec- 
ular dynamics simulations. In the 'glass' phase with no sym- 
bols, the relaxation time has become too large to be numeri- 
cally determined, the system is effectively a nonergodic glass. 
Note that the jamming transition at T — and ip = ip c can- 
not be accessed at thermal equilibrium. Fluids do not jam, 
only glasses do. 

is therefore also particularly interesting in the context 
of the question asked in our title, since the physics be- 
hind the density anomaly is a competition between en- 
ergy and entropy to minimize the free energy at thermal 
equilibrium. Briefly, the anomaly stems from the increas- 
ing difficulty, and thus the increasing entropic cost, to 
find states with little overlap, and thus with little en- 
ergy, when tp becomes large. Above tp*, the system thus 
prefers instead paying some finite energy to allow parti- 
cle overlaps, which can be done in many different ways 
and thus makes the structure of the fluid less ordered. 

As shown in Fig. 1, the jamming transition corresponds 
to an extreme case of the crossover observed at finite 
temperature, namely a transition betwen a law-tp phase 
were numerous states with no overlap can be found, and 
a laxge-tp phase where states without overlap only exist 
with a vanishing probability. This suggests theories of 
the liquid state fail to describe a sharp jamming transi- 
tion because they do not accurately describe this com- 
petition at very low temperatures and large densities, as 
we confirm in the following section. 



B. Numerical simulations 

Why does liquid state theory fail to account for the 
jamming transition? A direct answer is provided by nu- 
merical studies of the finite temperature behaviour of the 
model (1). In Refs. [35, 36] the dynamics of harmonic 
spheres was studied using molecular dynamics simula- 
tions. In Fig. 3 we report in the (tp, T) plane the state 
points for which thermal equilibrium could be reached, 
and the typical relaxation time measured with sufficient 
precision during the course of the numerical study. Nu- 



merical simulations thus indicate that it becomes increas- 
ingly difficult to reach equilibrium in the regime where 
density is large and temperature low, see Fig. 3. In the 
'glass' phase, structural relaxation does not occur dur- 
ing the timescale allowed by the numerical experiment, 
and the system is essentially frozen in a very long-lived 
metastable state. For practical purposes, it has all the 
characterisitcs of a glass [37], i.e. an amorphous (liquid- 
like) structure, which does not relax (solid-like) on the 
observation timescale. 

Therefore, simulations teach us that the free energy 
landscape of the system of harmonic spheres becomes 
very complicated in the glass regime [5]. It is this com- 
plexity which is responsible for the breakdown of liquid 
state theory, which cannot be used to explore the glass 
phase. 

In Fig. 3 we also note that along the T = axis, equi- 
libration cannot easily be achieved above tp sa 0.60, while 
most numerical determinations of tp c are much above this 
value, in the range ip c 0.64 — 0.66. Thus, we con- 
clude that even if one introduces thermal fluctuations 
into the game, the jamming transition cannot be crossed 
at thermal equilibrium because the glass transition inter- 
venes first, whatever the path in the (tp,T) plane which 
is followed, compressions or quenches. Therefore, the 
jamming transition can only be observed by compressing 
glasses — not fluids. 

The unavoidable conclusion is that a theory of the glass 
state is needed, instead of a theory of the liquid state, as 
we describe in the next section. 



IV. STATISTICAL MECHANICS OF 
HARMONIC SPHERES: GLASS PHASE 

Although the theoretical literature of the glass transi- 
tion is vast [37], there are not very many microscopic 
quantitative approaches, i.e. capable of formulating 
quantitative predictions starting from the knowledge of 
the interaction between the particles, as we attempt to 
do here. The mode-coupling theory of the glass transi- 
tion and the replica approach to the glass phase are the 
two examples we discuss in this section. 



A. Mode-coupling theory 

The mode-coupling theory of the glass transition was 
developed in the mid-80 's [38]. It is built using the tools 
first developed to describe the structure and dynamics of 
liquids, and, in its initial formulation, uses the formalism 
of projection operators to derive a closed set of dynamical 
equations of motion for time correlation functions of su- 
percooled liquids. In its common implementation mode- 
coupling theory can thus be seen as a 'black box' which 
is fed by structural information on the fluid (the two- 
point static structure factor of density fluctuations), and 
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FIG. 4: Dynamical phase diagram for harmonic spheres de- 
rived from mode-coupling theory [39] confirms the propensity 
of harmonic spheres to undergo a glass transition at low tem- 
peratures. The reentrant glass line at large ip is a dynamical 
counterpart of the density anomaly shown in Fig. 2. 



Fig. 4 had not been observed in Rcfs. [35, 36], because 
the range of volume fractions studied numerically was 
too narrow, but it was recently observed for a system of 
Hertzian spheres in Ref. [40]. 

Now, although the mode-coupling equations can in 
principle be used to describe the (nonergodic) dynam- 
ics in the glass phase, it needs as input the static struc- 
ture of the fluid, which is the very problem we ascribed 
ourselves at the beginning of the paper. Thus, although 
this approach is useful in determining the existence and 
location of the glass phase in the system, it can not be 
employed to describe the thermodynamic behaviour and 
structural evolution deep into the glass phase. To do so, 
we must turn to the thermodynamic replica approach. 

B. Thermodynamic replica approach 

1. The replica method 



provides as an outcome the time dependence of density- 
density autocorrelation functions at any wavevector. 

This short description shows that mode-coupling the- 
ory is in principle able to describe the features revealed by 
the numerical results shown in Fig. 3. Indeed we present 
in Fig. 4 results from the analytical study of the phase 
diagram of harmonic spheres. Part of these results were 
presented in Ref. [39], where all details can be found. 
In this work, a detailed analysis of the dynamic scaling 
properties near the fluid-glass transition occurring near 
T = was presented. In Ref. [40] a similar theoretical 
analysis was performed for an interaction potential very 
close to the harmonic potential studied here, where the 
harmonic exponent 2 in Eq. (1) is changed to 5/2, the 
so-called Hertzian sphere potential, yielding results qual- 
itatively similar to those presented in Fig. 4. 

Importantly the results of the mode-coupling analy- 
sis shown in Fig. 4 confirm the propensity of the har- 
monic sphere system to form a glass at large density and 
low temperature. However, the mode-coupling results are 
quantitatively not very accurate, as is well-known after 
two decades of mode-coupling studies in the field of the 
glass transition [37, 38]. 

Although quantitative agreement is not found between 
theory and simulations, mode-coupling theory makes one 
striking prediction for the phase behaviour shown in 
Fig. 4. Theory predicts the existence of a temperature 
regime where the fluid of harmonic spheres becomes a 
solid glass upon compression at constant temperature. 
Strikingly, upon compressing the glass further, there ex- 
ists a second critical volume fractions above which the 
glass is melted and becomes a flowing fluid. Such a 
glass melting at high density is typically not observed 
in standard models of liquids, and is the direct result of 
the particle softness [40]. We in fact already provided 
an explanation for this behaviour when we discussed the 
physics of the density anomaly in Sec. Ill A and Fig. 2. 
The reentrant glass transition predicted theoretically in 



The replica approach to the glass transition [26-29] can 
be seen as a modern and more systematic implementation 
of the self-consistent phonon theory developed long ago 
by Wolynes and coworkers [41]. The goal is to properly 
evaluate the thermodynamics of the glass phase, by care- 
fully taking into account the existence and proliferation 
of long-lived metastable states which are indeed respon- 
sible for the breakdown of liquid state approaches at low 
temperatures, as described in Sec. III. 

To account for the effect of metastable states, the par- 
tition function is decomposed into the respective contri- 
butions of inter- and intra-statcs free energies, 

z = J ' dfe -NPU-TW)) j ( 6 ) 

such that M(f) = exp[AT£(/)] represents the number of 
metastable states with free energy /, and £(/) is called 
'complexity'. This decomposition can be seen as a gen- 
eralization of the one based on the concept of inherent 
structures developed earlier by Goldstein [42] and Still- 
inger and Weber [43] . 

The replica approach is a computational tool devised 
to derive analytically the complexity of the system as it 
approaches the glass transition, and the thermodynamic 
properties of the system deep in the glass phase [26, 29]. 
It is based on the introduction of an effective temper- 
ature T e ff = T/m conjugated to the free energy / in 
Eq. (6), in the same way as temperature is conjugated to 
energy in standard computations. If m is an integer, it 
can be interpreted as a number of replicas of the original 
system [26]. A careful analysis shows that one is then 
able to deduce the thermodynamics of the glass from the 
thermodynamics of a m-times replicated liquid, with m 
being analytically continued to non-integer values. To 
compute the properties of the glass, one must then study 
the properties of an effective liquid which is a 'mixture' 
of m copies of the original system [27, 29]. We refer to 
Refs. [19, 44] for extensive reviews of a technique whose 



advantages and shortcomings, successes and failures, are 
by now well-established. 

As briefly mentioned in the introduction, this step is 
of course crucial to understand the physics of the glass 
transition, but it is equally fundamental to properly de- 
scribe the jamming transition which lies deep into the 
glass phase. The numerical protocols devised to study 
the jamming transition indeed all rely on a dynamical 
exploration of the ground state properties of the system 
at T = found either using gradient descent methods 
(which indeed produce inherent structures [43] ) , or rapid 
compressions. Given the topology of the phase diagram 
in Fig. 3, both rapid quenches and fast compressions hit 
the glass transition at some point, and from that point 
these numerical procedures in fact follow the zero tem- 
perature properties of long-lived metastable states across 
the jamming transition. The thermodynamic replica 
method, centered around the properties of these states, 
appears therefore very well-suited to describe the jam- 
ming transition in soft repulsive systems. 

As discussed at length in the glass literature [37, 44], 
the decomposition in Eq. (6) relies on the existence of 
infinitely long-lived states, such that a thermodynamic 
calculation makes sense. This is a typical mean-field as- 
sumption, because in finite dimensions metastable states 
only have a finite lifetime (which of course becomes ex- 
tremely large in the glass phase), and this makes the 
thermodynamic replica approach prone to criticisms. 

Approaching the jamming transition, though, the life- 
time of the states indeed diverges, so that we expect 
the mean-field approximation to behave better upon ap- 
proaching the jamming point than it does around the 
glass transition [19]. We conclude therefore that thermo- 
dynamic replica calculations offer a promising theoretical 
framework to study the ground state properties of har- 
monic sphere glasses. 
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FIG. 5: Thermodynamic, Tk(<p), and dynamic, Td(ip), glass 
lines determined from solving the replicated HNC approxima- 
tion are in good qualitative agreement with both simulations, 
Fig. 3, and mode-coupling results, Fig. 4. 



dynamic glass transition marked by the appearance of 
metastable states (it is closely related, in principle, to 
the mode-coupling singularity discussed above). 

These new results are shown in Fig. 5. The curves 
have the expected shape and suggest, once again, the ex- 
istence of a glass phase at low temperatures and large 
density, as seen in computer simulations, see Fig. 3. In 
particular, Td(tp) has the same qualitative behavior of the 
mode-coupling transition temperature in Fig. 4, suggest- 
ing that the reentrant glass transition scenario should be 
a robust feature of this system [40]. Note also that the 
ratio (T<2 — Tk)/Tk is often (inversely) correlated to the 
kinetic fragility of glassforming systems. Then, one pre- 
diction of the replicated HNC approach is that this ratio 
decreases on increasing density above ip ~ 0.6, which is 
consistent with the numerical observation that the liquid 
becomes more fragile when <p increases [35, 36]. 



2. First attempt: Replicated HNC 

In practice, the thermodynamics of the replicated liq- 
uid has to be computed using some liquid theory ap- 
proximation. The simplest of these approximations is 
the replicated HNC theory developed in Ref. [27]. It 
leads to a set of coupled equations for the diagonal and 
off-diagonal replica correlation functions, whose solution 
allows in principle to determine the full phase diagram 
of the system. However, it has been checked in the case 
of hard spheres that this approximation is correct in the 
liquid phase, while it fails badly in the glass phase and 
in particular close to jamming [27]. The origin of this 
failure has been discussed in [19, section IV], where the 
replicated HNC equations are presented in full detail. 

Therefore, this approximation can only be used to de- 
termine the glass transition lines Tk(<p) and Td(tp) for 
harmonic spheres. Within mean-field theory, the for- 
mer represents the thermodynamic glass transition where 
the complexity vanishes, while the second represents the 



3. Second attempt: The small cage expansion 

A more successful approach was later suggested by 
Mezard and Parisi, based on the physical idea that deep 
into glassy phases when temperature is small, the m repli- 
cas of the effective system will be undergo small vibra- 
tions within each state, suggesting that a perturbativc 
expansion of the free energy in the small cage size, A, 
can be performed. 

With this procedure, the free-energy of the replicated 
system is mapped onto the free-energy of a non-replicated 
system at effective temperature T e g = T/m. The glass 
phase is found to be described by values of m that are 
smaller than 1, so that the replicated liquid is found, 
in the small cage expansion, to be equivalent to a non- 
replicated liquid at higher temperatures. As a result of 
this purely analytical game, one can finally relate the 
thermodynamics of a state point of the original system 
located inside the glass phase, to those of a state point 
in the liquid phase of the effective system, where liquid 
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FIG. 6: The ideal glass transition line, Tk(<p), determined 
using the first order small cage expansion of the replica 
approach. Note that no transition is found below Tk ~ 
1.4 • 10 -4 , and the T — > limit cannot be studied. 



state theory, such as HNC, can reliably be applied. This 
scheme has been successfully applied to Lennard- Jones 
glasses in [28, 29]. 

We have studied the system of harmonic spheres us- 
ing the small cage expansion at first order. From the 
replicated free-energy, one can deduce the locus of the 
Kauzmann thermodynamic glass transition, and all ther- 
modynamic properties of the glass. We present our re- 
sults for Tx(ip) in Fig. 6, which is again in qualitative 
agreements with other schemes described above. 

Unfortunately, the small cage expansion, while being 
well-defined in the high density regime, fails in the jam- 
ming regime that we are interested in. The reason for 
this failure is that the non-analyticity at r — 1 of the har- 
monic spheres (which is mandatory for them to behave 
like hard spheres when T = 0) makes the free energy ex- 
pansion ill-behaved in the limit T — > when density is 
too small. This can be seen in Fig. 6, where the glass 
transition line cannot be followed down to Tk = 0, and 
the transition abruptly disappears at some volume frac- 
tion, which is an artefact of the small cage expansion. 

For the learned readers, this can be easily seen by look- 
ing at the first order expression given in Ref. [29]. The 
first correction to the free energy at small A has the form 
-/3AF = -12AP(p(m - 1) / °° r 2 g(r)AV(r)dr. We know 
that g(r) ~ e^ v M as T 0. We then deduce that 
-/3AF ~ 36A^/^Tp^^. Since the zcro-th order term of 

— f3F has a finite limit (the entropy of hard spheres), we 
conclude that the expansion is not defined in the limit 
T — > 0, because of the Ay/]3 prefactor, and thus cannot 



be used to study the approach to the jamming transition 
from above. 

This failure also explains why an alternative free en- 
ergy expansion was recently developed to study hard 
spheres [19]. In this approach, a different small cage ex- 
pansion was performed in powers of \J~~A (instead of A in 
the Mezard-Parisi scheme), consistent with the observa- 
tion that the A-expansion is divergent. Using the \f~A 
expansion, it was possible to study the jamming transi- 
tion (pressure, structure) on the hard sphere side [19]. 
However, this method only applies to hard spheres and 
it cannot be used to explore the jamming transition of 
harmonic spheres. 

Technically, a full description of the jamming transi- 
tion thus requires an expansion scheme which is able to 
describe the crossover between the yfA and the Ay/ft ex- 
pansions of the free energy, valid on both sides of the 
jamming transition. Such an new approximation forms 
the core of a separate publication [8] . 



V. CONCLUSION 

In this proceeding paper, we have justified our on- 
going effort to attack the purely geometric problem of 
soft sphere packing in three dimensions using the tools of 
equilibrium statistical mechanics. This philosophy thus 
suggests to add a temperature axis to the phase dia- 
gram and study the statistical mechanics of soft repul- 
sive spheres, the jamming transition being obtained in 
the limit of T — > 0. Adding thermal fluctuations to a 
situation where hard constraints need to be satisfied is 
actually a common tool, for instance in combinatorial 
optimization problems [5, 45]. 

However, we also showed that adding temperature does 
not immediately solve the problem, since one quickly re- 
alizes that jamming does not occur in a fluid but in the 
glassy part of the phase diagram where thermal equili- 
bration is not easily reached. Thus, describing jamming 
requires the development of accurate analytical tools to 
describe the structure and thermodynamics of soft repul- 
sive glasses, which is a delicate task. 

We showed that although conceptually feasible, pre- 
viously published analytical schemes actually fail near 
the jamming singularity, suggesting that new theoretical 
developments are needed to succcsfully derive a fully mi- 
croscopic, quantitative theory of the jamming transition. 
We shall report elsewehere the results of our work in this 
direction [81. 
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